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We study statistics of first passage inside a cone in arbitrary spatial dimension. The probability 
that a diffusing particle avoids the cone boundary decays algebraically with time. The decay expo- 
nent depends on two variables: the opening angle of the cone and the spatial dimension. In four 
dimensions, we find an explicit expression for the exponent, and in general, we obtain it as a root of 
a transcendental equation involving associated Legendre functions. At large dimensions, the decay 
exponent depends on a single scaling variable, while roots of the parabolic cylinder function specify 
the scaling function. Consequently, the exponent is of order one only if the cone surface is very close 
to a plane. We also perform asymptotic analysis for extremely thin and extremely wide cones. 
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I. INTRODUCTION 

Random walks are widely used to model natural pro- 
cesses in physics, chemistry, and biology (J-H. In partic- 
ular, first-passage and persistence statistics [a, @ of mul- 
tiple random walks underlie reaction-diffusion processes 
, spin systems [§4l0| , and polymer dynamics [ll|, [l2| . 

First-passage processes involving multiple random 
walks are equivalent to diffusion in a restricted region 
of space. For example, the probability that N ordinary 
random walks do not meet is equivalent to the probability 
that a "compound" walk in N dimensions remains con- 
fined to the region x\ < x 2 < ■ ■ ■ < x/v- This probability 
decays as t~ N ( N - 1 )/ 4 m the long-time limit [13l - tl5| . 

When there are only two or three particles, the com- 
pound walk is, in many cases, confined to a wedge, 
formed by two intersecting planes. Moreover, the well- 
known properties of diffusion inside an absorbing wedge 
[j| explain the long-time kinetics [l3l [l6l - [l9| . In general, 
however, the absorbing boundary is defined by multiple 
intersecting planes in a high-dimensional space. Apart 
from a few special cases, diffusion subject to such com- 
plicated boundaries conditions remains an open problem 

Our goal is to use cones in high dimensions to ap- 
proximate the absorbing boundaries that underlie such 
first-passage processes. In this study, we obtain analytic 
results for the survival probability of a diffusing particle 
inside an absor bing cone in arbitrary dimension. In a 
follow-up study 23], we demonstrate that cones provide 
useful approximations to first-passage characteristics of 
multiple random walks [2J |. 

We consider a single particle that diffuses inside an 
unbounded cone with opening angle a in spatial dimen- 
sion d (Figure [T]). The central quantity in our study is 
the probability S(t) that the particle does not reach the 
cone boundary up to time t. Regardless of the starting 
position, this survival probability decays algebraically, 
S ~ t~P , in the long-time limit. 

First, we find the exponent j3 analytically by solving 
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FIG. 1: Illustration of a cone with opening angle a. The 
initial location of the particle is parametrized by the radial 
distance r and the polar angle 0. 



the Laplace equation inside the cone. In dimensions two 
and four, this exponent is an explicit function of the open- 
ing angle a, and in particular, /3 = (n — a) /2a when 
d = A. In general dimension, we find j3 as a root of a 
transcendental equation involving the associated Legen- 
dre functions. 

Second, we derive scaling properties of the exponent. 
Interestingly, the exponent f3 becomes a function of a 
single scaling variable in the large-<i limit. We obtain the 
scaling function as a root of the transcendental equation 



D 2/3(2/ ) = with y = (cosa)Vd 



(1) 



involving the parabolic cylinder function D v . The ex- 
ponent ft is of order one only in a small region around 
a = it/2. The width of this region shrinks as d^ 1 / 2 in 
the infinite dimension limit. The exponent diverges al- 



gebraically, P{y) 
tially small, j3(y) 



II 



as y — ¥ 00, and it is exponen- 
\J y 2 /8ir exp(— y 2 /2) when y —> —00. 



Thus, in the large-d limit, the exponent (3 is huge if the 
opening angle is acute, and conversely, it is tiny if the 
opening angle is obtuse. Strikingly, if we fix the open- 
ing angle a and take the limit d — > 00, there are three 



2 



distinct possibilities, 



lim f3 d (a) 

d— >CG 



oo a < 7r/2, 
1/2 a = tt/2, 
a > tt/2. 



(2) 



Of course, a cone with opening angle a = tt/2 is simply 
a plane, and hence, f3(a = tt/2) = 1/2 for all d. 

Third, we study the limiting cases of very thin and 
very wide cones. The exponent diverges algebraically, 



when the cone is extremely thin. When the 



cone is extremely wide, the exponent is exponentially 
small, $ ~ (tt - a) d ' 3 . 

The rest of this paper is organized as follows. In Sec- 
tion II, we write the diffusion equation that governs the 
survival probability, and show that finding the leading 
asymptotic behavior of the survival probability requires 
a solution to the Laplace equation [25l - t29| . We present 
the solutions to this Laplace equation in two and four di- 
mensions in Section III, and for an arbitrary dimension 
in Section IV. The bulk of the paper deals with asymp- 
totic analysis for very large dimensions. In particular, 
we derive scaling properties of the exponent and obtain 
the limiting behaviors of the scaling function (Section 
V). Asymptotic results for extremely thin and extremely 
wide cones are detailed in Sections VI and VII, respec- 
tively. We also obtain the first-passage time (Section 
VIII) and conclude with a discussion in Section IX. 



II. THE DIFFUSION EQUATION 

Consider a particle undergoing Brownian motion |30L 
l3lj inside an unbounded cone in spatial dimension d. The 
opening angle a, that is, the angle between the cone axis 
and its surface, fully specifies the cone (Figure [T]). The 
range of opening angles is < a < tt, and for a = tt/2, 
the cone surface is planar. Moreover, the exterior of the 
cone is itself a cone with opening angle tt — a. In two 
dimensions, the cone is a wedge, and in three dimensions, 
the cone is an ordinary circular cone. 

At time t = 0, the particle is released from a certain 
location inside the cone. Our goal is to determine the 
probability that the particle does not reach the cone sur- 
face up to time t. By symmetry, this survival probability, 
S = S(r, 9, t), depends on the initial distance to the apex 
r, and the initial angle with the cone axis 9. Using a 
spherical coordinate system where the origin is located 
at the cone apex and the z-axis is along the cone axis, 
the pair of parameters (r, 9) are simply the radial and the 
polar angle coordinates of the initial location (Figure [1} . 

The survival probability fully quantifies the first- 
passage process. For example, the probability that the 
particle first reaches the cone surface during the time in- 
terval (t, t + dt) equals [— dS(r, 9, t) /dt] x dt. In general, 
the survival probability satisfies the diffusion equation Q 



dS{r,9,t) 
dt 



DV 2 S{r,9,t), 



(3) 



where D is the diffusion constant. The initial condi- 
tion is S(r,9,t = 0) = 1, and the boundary condition 
isS(r,a,t) = S(0,d,t) = 0. 

We are primarily interested in the large-time kinetics. 
Based on the behavior in two and three dimensions Q , we 
expect that the survival probability decays algebraically, 



S{r,9,t) ~ $(r,0)r 



(4) 



as t — > oo. The exponent f3 = [3d (ot) depends on the 
spatial dimension d and the opening angle a. The de- 
pendence on the initial location enters only through the 
amplitude $. 

We now substitute the leading asymptotic behavior (|4| 
into the diffusion equation ([3]). Since the time derivative 
becomes negligible in the long-time limit, the amplitude 
<!> satisfies Laplace's equation, 



V 2 $(r,(9) = 0, 



(5) 



subject to the boundary conditions $(r, a) — $(0, 9) = 0. 
The survival probability is finite and positive everywhere 
inside the cone, and consequently, the amplitude must 
be finite and positive, < $(r, 6) < oo for all 9 < a. In 
addition, d$/d9\g= = 0, to avoid a cusp along the cone 
axis. 

We use dimensional analysis [32j to solve Eq. (JSJ) - The 
survival probability is dimensionless, and hence, Q im- 
plies [$] = T' 3 where T denotes dimension of time. The 
amplitude $ depends on three variables: the radial coor- 
dinate r with dimension of length, [r] — L, the diffusion 
coefficient D with [D] — L 2 /T, and the dimensionless an- 
gle 9. As the quantity (r 2 / D) 13 is the only combination 
of the variables r and D with dimension T' 3 , we seek a 
solution in the form 



$(r, I 



tf>(0). 



(G) 



The angular function ?/>(#) must be finite and positive 
everywhere inside the cone, < ip(9) < oo for 9 < a, 
and it vanishes on the cone surface, i)(a) = 0. 

We now write the Laplacian operator in ([5]) explicitly 
by using 



V 2 = 



i_ 

dr'' 



d-1 d 



1 



d . d _ 2 d 
~£h- + r 2 {s\n9) d - 2 89 [Sm ' 89' 



Next, we substitute ^ into the Laplace equation and 
find that the angular function tp = ip(9) satisfies the 
eigenvalue equation (25l-|29| 



(sin( 



\d-2 



dO 



(sin 9) 



d-2 



dip 

~d~9 



+ 2(3(2/3 + d-2)ij; = 0. (7) 



The boundary conditions are ip'(0) = and ip(a) = 0. 

From equations ((4]) and ([6]), the leading asymptotic 
behavior is 



.S'( r.ti.l)^ r(0) (£1 



-0 



(8) 
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as t — >• oo. In particular, the survival probability grows 
algebraically with distance, S ~ r 2/3 . We note that the 
problem of finding the leading asymptotic behavior of the 
survival probability reduces to an electrostatic problem 
[iHliil as the amplitude satisfies Laplace's equation (|5|). 

Our goal is to find the exponent j3 and the angu- 
lar function ip{6). We expect that the exponent is a 
monotonically decreasing function of the opening angle 
a. Also, in all dimensions, (3(a = tt/2) = 1/2 because 
a cone with opening angle a — tt/2 is a half-space and 
consequently, the first-passage probability is identical to 
that of a particle in the vicinity of a trap in one dimension 
i- 



IV. GENERAL DIMENSION 

In general dimension, we transform the eigenvalue 
equation ([JJ using the variable /i = cos 9. In terms of 
this variable, the function -0 = ipd(fJ-) satisfies 



d/j, 2 v ' *~ dp, 



(d + 2/3(2^ + d- 2)ip = 0. (13) 



The boundary condition is ^(cosa) = 0. We now use a 
second transformation, 



V>(m) 



l 



M 2 ) 5/2 *(m) with 5 



d-3 



III. DIMENSIONS TWO AND FOUR 

We first discuss two special cases where explicit so- 
lutions are possible. In dimension two, the second or- 
der differential equation ([7]) reads ipgg + (2/3) 2 ip = 0, and 
the two independent solutions are simply cos(2/3(9) and 
sin(2/3#). The boundary condition ij)'(0) = excludes 
the latter and therefore, 



MO) = cos(2/30). 



Henceforth, the subscript indicates the dimension. Since 
the linear equation ([7J specifies the angular function 
only up to an overall constant, we set the prefactor to 
one throughout this paper. The boundary condition 
"0(a) = and the requirement that the angular function 
must be positive everywhere inside the cone together im- 
ply 2/3a = 7r/2. Therefore, the exponent is [B[ 



aw = Ta 



(10) 



As expected, the exponent is a monotonically decreas- 
ing function of a, and fli(TT/2) = 1/2. The exponent is 
minimal, yet finite, ^(tt) = 1/4, for an absorbing needle 
[35j . Thus, in dimension two, a diffusing particle reaches 
a needle with certainty. 

In dimension four, the transformation 
ip4(0) = (sm9)^ 1 u(8) reduces the eigenvalue equation 
([7]) to ugg + (2/3 + l) 2 tt = 0. Now, the two independent 
solutions are sin[(2/3 + 1)0] and cos[(2/3 + 1)9], but the 
latter is forbidden because the function ip must be finite. 
Therefore, 



M0) 



sin [(2/3 + 1) 
sin0 



(11) 



The boundary condition ip{a) — and the requirement 
tp(0) > for < a together give (2/3 + l)a = tt. Hence, 
the exponent is an explicit function of the opening angle 
in dimension four as well, 



ft (a) 



(12) 



The exponent vanishes, ft^a) — > as a — > tt, so it is 
no longer guaranteed that a diffusing particle reaches a 
needle. 



Substituting this form into (|13[) shows that the auxil- 
iary function ^ = ^d(fJ>) satisfies the associated Legendre 
equation [36l. |37| 



{1 -^d^ 



2/i — 
d[i 

(2/3 + 8) (2/3 + 5 + 1) 



1-fi 2 



(14) 
* = 0. 



(9) 

The two independent solutions of this equation are 



P, 



(/i) and Q^p+si^), the associated Legendre func- 



2/3+6 

tions of degree 2/3 + 5 and order S. In even dimensions, 
the first solution is not physical because it implies a diver- 
gence of ijj as — > [37|. In odd dimensions, the second 
solution is excluded for the same reason. Therefore 13811 . 



J(sin0)- 5 P^ +5 (cos0) dodd, 

ipd{0) = i / . \-Sr>6 i a\ a 

[(sm0) Q2/3+s{ cos 0) a even. 



(15) 



The boundary condition ip(a) = relates the exponent 
/3 and the opening angle a, 



P. 



2/3+5 



(cos a) 







Q S 2p+s( cosa ) = 



d odd, 
d. even. 



(16) 



In general dimension, the exponent /3 is the smallest root 
of the transcendental equation (|16[) involving the asso- 
ciated Legendre functions. We must always choose the 
smallest root because ip(0) > for all < a. 

We can verify that the exponent /3 decreases monoton- 
ically with the opening angle a, and that (3(tt/2) = 1/2 
(Figure [2]). In all dimensions, the exponent diverges for 
extremely thin cones, a — ¥ 0, and when d > 3, the expo- 
nent vanishes when a — > tt. 

Let us fix the opening angle and increase the dimen- 
sion. There are two possibilities: (i) if a < tt/2, the 
first-passage process speeds up with increasing dimension 
because (3 increases (Figure [5]) and ip(0) declines (Figure 
[3k); (ii) if a > tt/2, the first-passage process slows down 
because /3 shrinks (Figure [2]) and ip(0) grows (Figure [3|d). 
Ultimately, in the infinite-dimension limit, first passage 
becomes instantaneous for acute angles but infinitely long 
for obtuse angles. 
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FIG. 2: The survival exponent j3 versus cos a at different 
dimensions. This exponent is given explicitly by (|10p and 
(|12|l in dimensions d = 2 and d = 4, respectively. In general, 
the exponent is a root of the transcendental equation (|16|1 . 



In dimension three, ^ = 0, and the solution (|15p re- 
duces to the Legendre function of index 2/3, namely, 
■03(0) = P2p(cos8) [H, H3- Also, the exponent /3 is the 
smallest root of the transcendental equation [5j] 



P2fi(cosa) = 0. 



(17) 



For half-integer values of the exponent, j3 = n/2, the 
angular function ipd(9) is a polynomial of degree n in 
cos 8. This follows directly from equation (fl~3|) . For ex- 
ample, 



cos 
= <( dcos 2 6»- 1 
(d + 2) cos 3 9 



3 cos 



= 1/2, 

= 1, 
= 3/2. 



(18) 



In three dimensions, the polynomials coincide with the 
Legendre polynomials: ip3{0) = P n (cos9) when (3 = n/2. 
Using i)(a) = 0, we find the opening angles for which the 
exponent is a half-integer: 





f;/2 


cos a 


Pd (a) = < 




cos a 




^3/2 


cos a 



o, 

l /Vd, 
^3/(d + 2) 



(19) 



V. SCALING PROPERTIES 



The special values listed in p9|) suggest that the sur- 
vival exponent has the scaling form 



Pd{ct) — > P{y) with y — (cosa)Vd, 



(20) 



in the limit d — > oo. From equation (|19p . we deduce 
/3(y = 0) = 1/2, /3(y = 1) - 1, and /% = V§) = 3/2. To 
show the scaling behavior (|20[) in general, we introduce 
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FIG. 3: The function ip(9) versus the normalized polar angle 
9/a at different dimensions. The function ip(6) is normalized 
such that ^(u) = 1. Figure 3a shows the behavior for an acute 
opening angle, a — n/4, and Figure 3b shows the behavior 
for an obtuse opening angle, a — 3n/4. 



the variable z = u Vd. Performing this scaling trans- 
formation on the Laplace equation (fT3")) shows that the 
angular function depends on a single scaling variable, 
ipd{^) ~~ V'( z )i m the d — ¥ oo limit. The scaling func- 
tion ip{z) satisfies 

+ 20^ = 0. (21) 

The boundary condition is ip(z = y) = 0, and addition- 
ally, tp(z) > for all y < z < oo. 

Using the transformation ip{ z ) = exp(z 2 /4)zt(z), we 
recast (|2"T]) into the parabolic cylinder equation [39| 



2/3- 



u = 0. 



The two independent solutions of this equation are 
D2p(z) and D2p(—z) where D v (z) is the parabolic cylin- 
der function of index v. The second solution is not 
physical: the asymptotic behavior D v (—z) ~ exp(z 2 /4) 
as z — > oo 39] implies a divergent survival probability in 
the limit 9 — > 0. Therefore, 



1>{z) = e z2 ^D 2 p(z). 



(22) 





FIG. 4: The exponent /?, specified by equation (|23p . versus 
the scaling variable y = (cosa)Vo. 



As always, we set the overall proportionality constant to 
one. The boundary condition ip(z = y) = relates the 
scaling function /3(y) and the scaling variable y, defined 
in (EOD, 



(23) 



The proper solution is the largest root of the parabolic 
cylinder function. For half-integer values of the expo- 
nent, the parabolic cylinder function is related to a Hcr- 
mite polynomial, D 2 /3(z) cx H n (z/V2) cxp(— z 2 /4), and 
hence, (f23jl is equivalent to H n (y/V2) = 0, where the 
largest root is the appropriate one. Therefore, the half- 
integer values of the scaling function /3(y) occur at ze- 
roes of Hcrmitc polynomials. In addition to the afore- 
mentioned values, we also quote /3(v 3 + \/6) = 2 and 

P(y/5 + Vw) = 5/2. 

The scaling behavior (|2T))) is valid in the limits a —} tt/2 
and d — Y oo with the scaling variable y — (cosa)y/~d or 
alternatively y = (tt/2 — a)Vd kept finite. Hence, the 
scaling function in (|2"31 quantifies the shape of (3 in a 
"scaling window" [401442] centered ona = tt/2. The size 
of this window shrinks as <i -1 / 2 at large dimensions. As 
shown in Figure^! /3 vanishes as y — » — oo, and it diverges 
as y — > oo. Surprisingly, if we fix the opening angle a 
and then take the d — > oo limit, there are three distinct 
possibilities, as stated in @: the exponent vanishes if 
a > 7r/2, it always equals 1/2 when a = tt/2, and it 
diverges if a < tt/2. 

We also note that the approach to the scaling behavior 
is not uniform. The scaling exponent /? converges rapidly 
to the scaling function for negative y, but the convergence 
is quite slow for positive y (Figure [5]) . 

To find the asymptotic behavior when y — > oo, we use 
the fact that the largest root £ = £(z/) of the parabolic 
cylinder function of index v, D u (£) = 0, is located at [36[ 



(4^-1^-1/6, 



FIG. 5: The convergence to the scaling behavior. The expo- 
nent /3 = Pd (a), specified by (|16[) . is plotted versus the scaling 
variable y = (cosa)\/d for d = 10 and d — 100. Also shown 
is the scaling function /3(y). 



when v — > oo. Here, a\ = —2.338107 is the first root 
of the Airy function. Substituting v — 2(3 into this ex- 
pression, we find y ~ (8/3) 1 / 2 - |ai | (2/3) -1 / 6 as j3 -> oo. 
Therefore, the leading behavior is j3 ~ y 2 /8. Further- 
more, the first correction to this asymptotic form is 



P~ ^- + \ ai \2- 5 / 3 



2/3 



(24) 



Thus, the exponent j3 diverges algebraically when 
y — > oo. 

We use perturbation analysis to find how j3 vanishes in 
the complementary limit y — > oo. Equation (|2ip shows 
that the angular function becomes constant when /3 — > 0. 
Thus, ip(z) ~ l + f3g(z). Substituting this form into ([2T]). 
the correction function g = g{z) obeys 



9zz - zg z + 2 = 0. 
The boundary condition 

5 ,(0) = V2t7 



(25) 



(26) 



follows from the small-/?, small- z, behavior of the scaling 
solution 1221 



i){z) ~ 1 + 



= 1- 



v/2r(i-^) 
r(-/3) 



Here, we used the values D 2 p(0) = V?r2^/r(l/2 - 0) 
and L>^(0) = -V^2' 3 + 1 / 2 /r(-/3) j^, as well as the 
identities T(l/2) = ^ and (-/3)r(-/3) = T(l - /3). 

Using the integrating factor I(z) = exp(— z 2 /4), we 
simplify Eq. (|23|) to d(g z I)/dz — —21. Integration of this 
equation subject to the boundary condition (|26p gives 



j(z) = 2 dte t2/2 

Jo Jt 



ds t 



-s 2 /2 



(27) 
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The second integral appr oaches a constant, and 
consequently, g{z) ~ — ^J&k/z? exp(z 2 /2), when 
z — > — oo. We now impose the boundary condition 
ip(y) = 1 + g{y) = 0, and find that the exponent is 
exponentially small, 



(28) 



in the limit y — > — oo. 

In summary, the scaling function j3(y) has the following 
extremal behaviors 




\/y 2 /S>n exp(-y 2 /2) y -> -oo, 
y 2 /8 y^oo. 



(29) 



d 


B d 


2 


0.785398 


3 


1.202412 


4 


1.570796 


5 


1.915852 


6 


2.246704 


7 


2.567811 


8 


2.881729 


9 


3.190080 


10 


3.493966 



TABLE I: The coefficient B d in (32) for d < 10. Equations 
(9]) and (JTTJ) yield f?2 = 7r/4 and B4 = tt/2, respectively. 



At large dimensions, these limiting behaviors apply only 
inside the scaling window, that is, when |tt/2 — a\ is of 
the order d~ 1 ^ 2 . In the next two sections, we perform 
asymptotic analysis for the limiting cases of extremely 
thin cones (9 — > 0) and extremely wide cones (9 — > tt). 



VI. THIN CONES 



The explicit expressions (|T0|) and f| 12[) show that the 
exponent is inversely proportional to the opening angle, 



a 



when the cone is extremely thin. We anticipate 



that this divergence is generic. 

When the cone is very thin, we have sin 9 
tion @ simplifies, 



d 2 i> 



d— 2 dib , nm , , 

d 9 2 i + w^=°- 



9 and equa- 



(30) 



This equation holds as long as da 2 C 1. In writing this 
equation we tacitly assumed a divergent /3. We now intro- 
duce the scaling variable x — 2/39 and transform equation 
(1501) as follows, 



^Pxx 



d-2 



■ i>x + i> = 0. 



Next, we seek a solution in the form ip(x) 



where as in (fT5|) , 6 = ■ Performing this second trans- 
formation gives the Bessel equation [43( 

x 2 u xx + x u x + (x 2 — <5 2 ) u = 0. 

This equation has two independent solutions: the Bessel 
functions Js(x) and Yg{x). The boundary condition 
ip'(0) = excludes the latter and hence, 



ip(x) = x J$(x). 



(31) 



The other boundary condition, ip(a) = 0, gives an im- 
plicit relation Js(2fia) — between the exponent /3 and 
the opening angle a. We see that the exponent is in- 
versely proportional to the opening angle in all dimen- 
sions: 



(i d {a) ~ B d a~ 



with Bri 



has), 



(32) 



as a — >• 0. Here, C(<5) is the first positive zero of the 
Bessel function, Js{C) = 0. Table I lists the coefficients 
B d for d < 10. 

At large dimensions, we use the asymptotic behavior 

aS)o,S+\a 1 \(S/2) 1 / 3 

when 6 — > 00. Again, a\ = —2.338107 is the first root of 
the Airy function. Thus, to leading order, the prefactor 
Bd is linear in the dimension, Bd — d/A. The correction 
to the leading asymptotic behavior is given by 

B d ^^ + |a 1 |2- 5 / 3 d 1 /3 j (33) 



for d — > 00. The divergence (|32|) shows that a diffusing 
particle quickly reaches the cone surface when the cone 
is thin. Moreover, equation (|33p implies ft ~ d/Aa, and 
hence, the first-passage process becomes even faster as 
the dimension increases. 



VII. WIDE CONES 

For all d > 3, the solutions to equation (flrj|) show that 
the exponent vanishes when the interior of the cone oc- 
cupies all of space (Figure [2]). In this case, equation (7]) 
becomes 



1 



(sin(9) d - 2 d6 



' d9 



+ 2/3(d - 2)ip = 0. (34) 



We obtain j3 by repeating the perturbation analysis lead- 
ing to (2SJ. For all d > 3, we write ij)(6) ~ 1+0 g{6). The 
boundary condition ip'(0) = implies g'(0) = 0. From 
(l34t . the correction g = g(9) obeys 



d 



(sin( 



' d9 



-2(d- 2) (sin l 



\d-2 



(35) 



Integrating this equation twice and using the boundary 
condition g'(0) = 0, we find the correction up to a con- 
stant, 

g(9) = g - 2(d-2) £-J*— f\ (sin^ 2 . (36) 
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d 


3 


4 


5 


6 


7 


8 


9 


10 


A d 





i 

2tt 


i 

4 


l 

7T 


3 
8 


4 

3-jr 


15 
32 


8 

5tt 



TABLE II: The coefficient A d given by equation for 
d < 10. 



In the limit 9 — > it, the second integral approaches the 
constant 



, 2 r(^i)r(|) 



S d = / d6 (sin t 



r(|) 



Moreover, the term (sin</>) _ ( d ~ 2 ) ~ (7r — cf>)~( d ~ 2 ' > domi- 
nates the first integral and therefore, the correction func- 
tion diverges, 

9(o)*- 2{d d Z 2 l Sd (*-er ld - 3) > 

as 9 — > it. We now impose the boundary condition 
ip(a) = 1 + (3g(a) = and find that the exponent /3 
vanishes exponentially as a — > n, 



/3 d (a)~A d (Tt-a) d - 3 . 

1 — d — 3 
ld - 2(d-2)S d 



The coefficient A d = ^,,t is given by 



A d 



r(¥) 



2T 



(37) 



(38) 



In particular A^ = 1/2tt, in agreement with equation 
(fTU). Table II lists the coefficients A d for d < 10. The 
asymptotic property T(x + a)/T(x) — > x a as x — > oo 
shows that the coeffic ient grows algebraically with di- 
mension, Ad — v d/8n. 

In the marginal case <i = 3, the coefficient in (|3"T|) van- 
ishes, A 3 = 0. In this case, the exponent vanishes gently, 



4 In 



2 

7T— a 



(39) 



as follows from ([TT]) 43[ . This behavior is consistent with 
the near-7r behavior of (|36|) . g(0) ~ — 41n^g. 

As the opening angle approaches its maximal value, the 
cone surface turns into an infinitely long, yet infinitesi- 
mally thin needle. In two dimensions, a diffusing particle 
is bound to reach such a needle. Yet, in dimensions three 
and higher, the particle may or may not reach the needle. 
Moreover, the exponentially small exponent <f3T|) shows 
that the first-passage process becomes extremely slow at 
high dimensions. 



VIII. FIRST-PASSAGE TIME 

We now briefly discuss the first-passage time. Let 
T(r, 9) be the average duration of the first-passage pro- 
cess, namely, the average time it takes a particle released 



at (r, 9) to reach the cone surface for the first time(Fig- 
ure[T]). This quantity obeys the Poisson equation 



DV 2 T(r,/ 



1, 



(40) 



and the boundary condition T(r,a) = T(0,9) = 0. 
As in (|6]), we use dimensional analysis and write 
T(r,0) = (r 2 / 'D)U(9) where U(9) is a dimensionless 
function of the angle. From Eq. (j40|) . we get 



1 



(sin6») d - 2 dO 



(sm9) d - 



,dU 

~dJ 



2dU 



-1, 



(41) 



while the boundary condition becomes U (a) = 0. 
The linear equation (|41|) has the particular solution 
U = —l/2d, and thus, we seek a solution in the form 
U{9) = u{9) - l/2d. The function u = u{9) obeys the 
homogeneous linear differential equation 



1 



(sin6») d - 2 dO 



{sm9) d - 



, du 

~d~9 



2du = 0, 



subject to the boundary condition u(a) = l/2d. This 
equation is a special case of (O, and using ([LB]) , we im- 
mediately find u(9) — C(d cos 2 9—1) where C is set by 
the boundary condition. Finally, we find that the first- 
passage time has the compact form 



TM = ^ C0S29 



cos 2 a 



(42) 



2D dcos 2 a-l 
This first-passage time is finite if and only if fj > 1. 

IX. DISCUSSION 

In summary, we studied first-passage kinetics for a par- 
ticle diffusing in a cone. In all dimensions, the probabil- 
ity that the particle does not reach the cone boundary 
decays algebraically with time. We found the exponent 
underlying this power-law behavior as a root of a tran- 
scendental equation involving associated Legendre func- 
tions. We also obtained scaling and extremal properties 
of the exponent. 

Our results generalize the known properties of first pas- 
sage in two- and three-dimensional cones [5[ to arbitrary 
dimensions. Moreover, the statistical physics perspec- 
tive, where scaling plays a central role, extends rigorous 
studies of diffusion in cones in probability theory [25l — 
l28j and potential theory [2!|. Scaling implies that the 
behavior in a narrow window becomes universal: cones 
with different combinations of a and d, have the same ex- 
ponent /3, as long as the scaling variable y = (cos a)Vd 
is the same. The exponent is exponentially small when 
y — > -co, and it grows algebraically when y — > oo. 

As the dimension increases, the first-passage process 
speeds up if the opening angle is acute but slows down if 
the opening angle is obtuse. This behavior is reflected by 
the asymptotic behavior for fixed 9 in the large-d limit. 
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We merely quote the results of the corresponding asymp- 
totic analysis, 



Pd(a) 



f X 
\ sin a 



i 



C(sincv) d 



a < 7r/2, 
a > tt/2. 



(43) 



The exponent is proportional to the dimension at acute 
angles, and it decays exponentially with dimension at 
obtuse angles. The asymptotic behavior (|43l) is consistent 
with the limiting behavior of the scaling function (|29|) as 
well as the asymptotic results (f3"2"j) and (|3"T|) . 

Kinetics of first-passage in a cone provide useful in- 
formation about first-passage problems involving mul- 
tiple random walks. For example, let us consider the 
random-walk problem mentioned in the introduction. 
The probability that the trajectories of N ordinary ran- 
dom walks do not meet up to time t is equivalent to 
the probability that a compound random walk in N di- 
mensions remains inside the region x\ < x% < ■ ■ ■ < xm- 
Since there are N\ permutations of the positions, the 
total solid angle inside this "allowed region" accounts 
for a fraction 1/N\ of space. We replace the allowed 
region with an A^-dimensional cone that has the same 
solid angle, a N_1 ~ 1/N\. Using the Stirling formula, 
Nl ~ y/2nN(N/e) N , the opening angle is a ~ e/N. The 
survival probability decays algebraically, S(t) ~ t~^ N , 
and using (|52 ]) -(|3"3" 1) . the exponent is fi N ~ N 2 /ie. This 
cone approximation correctly gives the asymptotic N- 



dependence as the exact exponent is N(N— 1)/4 [13l - ll5l |. 

First passage in the exterior of this narrow cone gives 
an approximation for the probability that the order of 
N random walks does not turn into the mirror image 
of the initial state. If the initial ordering of the ran- 
dom walks is {1, 2, • • • , N - 1, N}, then S(t) is the prob- 
ability that the particles do not reach the configuration 
{N, N — 1, . . . , 2, 1} up to time t. Since the angle of the 
cone is now n — a ~ e/N, the asymptotic behavior (|37[) 
gives S(t) ~ t~P N with the tiny exponent 



/3 N ~ y/N/8n (e/N) 



N-3 



In a follow-up study, we use cones to understand other 
random walk problems (23j. 

With a few exceptions such as paraboloids [3, [45| , the 
study of first passage in elementary geometries is still in 
its infancy. We considered first passage in basic circular 
cones. A natural extension is to generalized cones [46J, 
defined as domains with the property that all rays ema- 
nating from the apex do not intersect the cone boundary. 
Understanding first-passage properties in such general- 
ized cones ultimately requires a solution to the Laplace- 
Dirichlet boundary value problem. 

We thank Sidney Redner for useful discussions. This 
research has been supported by DOE grant DE-AC52- 
06NA25396 and NSF grant CCF-0829541. 
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